Laboratory measurements of electrostatic solitary structures 
generated by electron beam injection* 
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Electrostatic solitary structures are generated by injection of a suprathermal electron beam par- 
allel to the magnetic field in a laboratory plasma. Electric microprobes with tips smaller than 
the Debye length (Azj e ) enabled the measurement of positive potential pulses with half-widths 4 to 
25A.D e and velocities 1 to 3 times the background electron thermal speed. Nonlinear wave packets of 
similar velocities and scales are also observed, indicating that the two descend from the same mode 
which is consistent with the electrostatic whistler mode and result from an instability likely to be 
driven by field-aligned currents. 

PACS numbers: 52.35.Sb, 52.35.Mw, 52.72,+v, 94.05. Fg, 52.35.Fp 
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Electrostatic solitary structures identified as 
Bernstein- Greene-Kruskal electron holes p] have 
been abundantly observed in various active regions of 
space subject to large-scale current systems or energy 
dissipation (@, S| and references therein). This includes 
regions where magnetic reconnection occurs [J, [5[ . Elec- 
tron holes were also recently detected in a laboratory 
magnetic reconnection experiment @. Electron holes 
can be generated by energetic streaming electrons and 
are thought to play an important role in scattering these 
electrons. However, although these small-scale (one to 
tens of Debye lengths, Arj e ) solitary structures seem 
ubiquitous in key regions of space their exact origin 
often remains unclear. In the laboratory, experiments 
dedicated to electron holes were carried out in a strongly 
magnetized Q-machine Q- These holes were generated 
by a voltage pulse and had sizes comparable to the 
plasma column radius, making comparison with holes 
observed in space difficult. 

This Letter reports measurements of electrostatic soli- 
tary waves generated by an electron beam injected into 
a magnetized low-/? plasma column much larger than the 
structure scales (Fig. [T]). The experiment was conducted 
at the upgraded Large Plasma Device (LAPD) Q at the 
University of California, Los Angeles. The helium plasma 
column has a 60 cm diameter, is 17.1 m long and pulsed 
at 1 Hz with pulses lasting several milliseconds (Fig. [TJt) . 
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TABLE I: Main plasma parameter range in the afterglow 
phase during which the measurements are taken. Magnetic 
field strength and beam energy are fixed parameters for each 
experiment, while densities and temperatures change with 
time in a highly reproducible way. 

Magnetic field strength, B 100-750 G 

Background electron density, n e l.l-5.7xl0 9 cm -3 

Background electron temperature, T e 0.18-0.20 eV 
Electron beam voltage 60-120 V 

Debye length, Au e 44-95 /im 

Electron mean-free path, A ee 0.11-0.42 m 

Thermal electron gyroradius, r ge 13-106 /an 
Plasma frequency, f pe 297-681 MHz 

Electron gyrofrequency, f ce 0.53-7/ pe 



An electron beam 0.4 to 1 cm in diameter is injected from 
a 3 mm diameter LaB6 crystal source for about 140 fis 
along the axis of the column in the afterglow phase, be- 
tween 50 and 150 ms after the end of the discharge pulse. 
The beam density 5 cm from the source is approximately 
25% of the background electron density. The magnetic 
field strength, plasma density and beam voltage can be 
changed from experiment to experiment. The range of 
the main plasma parameters is summarized in table |U 
Floating potentials and electric fields are measured at 20 
GHz (38fp e for the average plasma density) by a micro- 
probe Q with 10 /im wide tips (at most Ad c /4) and 
with separations ranging from 40 to 130 /im (~ Xdc) 
[Fig. QJd]. Microprobe tips smaller than Xu e are critical 
to study Debye-scale structures without destroying them. 
The probe is located on the column central axis, 6-24 cm 
or 150-5600Az) e away from the beam source and parallel 
to the magnetic field and the column axis. 

The fluctuations excited by the beam typically have 
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FIG. 1: a) Schematic of probe and small electron beam in the 
LAPD device, b) Highly magnified photograph of the electric 
field MEMS microprobe. The conductor is made from gold 
and the insulating material is polyimide. The dimensions are 
w=10 fim, g = 30 fim and d= 70 /im. The wires are coaxial 
and tapered in their expansion regions to maintain a 50 f2 
impedance 



a broadband spectrum ranging from the bottom of the 
amplifier's frequency range (20 MHz) to the electron 
plasma frequency f pe and above. For the example shown 
in Fig. the power-spectrum appears exponential for 
0.2 < f / f pe < 0.8. Fluctuations in this range display a 
high degree of phase coherence between signals measured 
on probes separated by 2.3A_o e (Fig. [5b), and consist of 
a mixture of quasiperiodic waves, wave packets and soli- 
tary pulses. Solitary waves alone are also known in some 
cases to result in exponential spectra [loj . 

An example of a solitary positive potential pulse is 
shown in Figs. [2pde. The figure displays the potentials 
recorded on two probe tips separated by A = 1.77 Xdc 
(100 /im) along the magnetic field direction, the elec- 
tric field inferred as E = —Sep/ A and a lower bound 
(in magnitude) to the actual density perturbation cal- 
culated as p — eodE/dx assuming a slablike structure 
propagating at a constant velocity. Quantities are made 
dimensionless using the plasma parameters in the un- 
perturbed state before beam injection. In this example 
n e = 3.5 x 10 9 cm- 3 , T e = 0.2 eV and B = 750 G 
(fpe/ fee — 1/4). The pulse conserves its shape while 
traveling between the probe tips. The perturbation is 
small in the sense that ecj> < T e and An e <C n e . The time 
delay between probes is 0.2 f~^ (0.3 ns), implying a par- 
allel velocity in the same direction as the beam electrons 
W|| = 1.98w Te = 371 km/s with v Te = (T e /m e ) 1/2 , and 
a parallel half-width L\\ — 10.3A_De (575 /im). The pulse 
has a symmetric positive potential with dipolar electric 
field and a charge density which is positive at the center 
and negative near the edges. These properties are consis- 
tent with a Bernstein-Greene- Kruskal electron hole [lj. 




FIG. 2: Power spectrum (a) and coherence (b) for potential 
time-series on two field-aligned probes separated by 2.3Aoe. 
A solitary potential pulse (c), its bipolar electric field (d), 
and estimated charge density perturbation (e). In this and 
all other figures the parameters used for normalization corre- 
spond to the unperturbed state before beam injection. 
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FIG. 3: Scatter plots of solitary pulse half-width versus am- 
plitude. The shaded area corresponds to a theoretically inac- 
cessible region frD ]. Points are color coded according to the 
magnetic field ranging from 100 G (black/dark red) to 750 
G (light orange). The dotted lines indicate the median val- 
ues. The inset shows a probability histogram of the structure 
velocities. 
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An automated search for isolated peaks in the electric 
potential time series for plasma parameters in Table U 
yielded a total of 363 solitary structures over 35 experi- 
ments with various beam voltages, plasma densities and 
magnetic fields. The propagation time from one probe 
tip to another is determined by cross-correlation analy- 
sis, and velocities cross-checked between different pairs of 
probe tips. All structures propagate in the same direction 
as the beam and display a main positive potential peak. 
40% have symmetric dipolar electric fields, while 30% are 
more tripolar and the remainder mostly dipolar with a 
notable asymmetry. Scatter plots of their main parame- 
ters are shown in Fig. [3] The half-widths range from 4.5 
to 24.4A.De with a median of 10.2Ar>e- Amplitudes are all 
smaller than T e /e with a median of 0.13T e /e. 80% of the 
structures have velocities between 1.3 and 2.3vre (Fig. 
[3j inset). 

The typical scales are comparable to those usually re- 
ported in space observations, albeit a bit larger. Ref- 
erences @, [|[ report average half-width slightly larger 
than A.D e and velocity slightly smaller than VTe- This dif- 
ference might be partly explained by the normalization, 
which in our cases uses background parameters measured 
prior to the beam injection. In particular the plasma in 
which the solitary waves propagate is likely warmer due 
to the interaction with the beam, making the structure 
velocities closer to the actual thermal speed. Indeed it 
was predicted that electron holes faster than 2vr e are 
unstable [l( and are therefore unlikely to be observed. 
Normalization questions apart, electron hole sizes signif- 
icantly larger than A^e are theoretically allowed, elec- 
tron holes of parallel size 60X]j e have been reported in a 
laboratory experiment @. 

The widths and amplitudes do not display any one- 
to-one relation (Fig. [3]), unlike Korteweg de Vries-type 
solitons. However the positivity of the trapped electron 
distribution imposes an upper limit to the hole ampli- 
tude for a given scale ll|. All the amplitudes lie well 
within the allowed region for stationary Id (infinite per- 
pendicular size L±) or magnetized 3d (with finite L±) 
electron holes predicted by The electron trapping 

frequency f tr = (e/c|| \E\\ \/m e ) 1 t 2 is found to be typically 
a few percent of / ce , showing that electrons remain mag- 
netized within these structures. However parallel sizes or 
the amplitudes have no clear relationship to the magnetic 
field intensity (color-coded in Fig. [3]). 

Besides electron holes, irregular fluctuations and wave 
packets are observed. Fig. H^b shows an example of a 
solitary structure shortly followed by a wave packet. The 
approximately constant phase-shift shows that the soli- 
tary wave and wave packet have very similar velocities 
and scales and therefore might be related to the same 
plasma mode. Measurements on different probes allow 
estimation of the local wave-number and spectral den- 
sity based on the instantaneous cross-spectra The 
density and the average experimental dispersion relation 
w(/c||) are shown in Fig. The fluctuations follow a 
weakly dispersive curve, although the spectral density 
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FIG. 4: (a) potentials on two probes showing a solitary struc- 
ture, (b) a wave packet shortly following the previous solitary 
structure, (c) experimental dispersion relation oj(k^) (the 
solid blue curve is the average), (d) phase velocity derived 
from the experimental dispersion relation, (e) 8kB derived 
from Eq. [T] The horizontal dotted curve in (c) is the in- 
verse median time-scale of the electron holes in this particu- 
lar experiment, and the one in (d) is their median velocity. 
The vertical dotted line in (c)-(e) is the wave- number corre- 
sponding to their median parallel width. In this experiment 

/ce/ fpe — 3. 



becomes smaller and more spread out as the frequency 
increases toward f pe . The average experimental disper- 
sion curve extends to the temporal and spatial scales 
where the electron holes are observed and implies par- 
allel phase (and group) velocities about 2vTei very simi- 
lar to the measured velocities of the holes (their median 
velocity for these experimental parameters is 2.lvxe, the 
horizontal dotted line in Fig 4d). In this frequency and 
wave-number range the dispersion is consistent with the 
whistler mode at the resonance angle, which is essentially 
electrostatic and has a lower cutoff at the lower-hybrid 
frequency u>lh ■ In the frequency and wave- number range 
we observe the waves ion and beam contributions can be 
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neglected and their dispersion relation reduces to 



where P and S are standard elements of the cold plasma 
dielectric tensor (l3j . Electrostatic whistler mode waves 
following the average experimental curve w(fcy) have an 
angle predicted by Eq. [T] shown in Fig. [4J5. The waves 
are predicted to become more perpendicular as k\\ de- 
creases and have k±Xne between 0.33 and 0.39. Electron 
holes and electrostatic whistler mode waves known as 
VLF saucers have been found to be closely associated in 
the auroral regions [T3|. They have also been observed 
together in the polar cusp [HI, [Hj]. The close associa- 
tion between the electron holes and electrostatic whistler 
mode waves suggests that both originate from the same 
source, or that one produces the other [l7| . 

Knowledge of the electron distribution function is key 
in sorting out the type of instability that is responsible 
for the generation of the observed waves and perhaps the 
solitary structures. The beam distribution was measured 
using a swept Langmuir probe capable of 100 V sweeps. 
The derivative of the I-V characteristic curve is propor- 
tional to the distribution function. 32 cm away along the 
background magnetic field from the 66 V beam source no 
beam is visible, and only a 35 eV tail is left. However at 
the probe location where most of the solitary structures 
are measured (about 5-6 cm from the beam source), a 
warm beam with a mean energy of 62 eV, thermal energy 
Tft w 5 eV and a density of approximately 25% of the to- 
tal density is observed. In addition a high energy tail on 
the background distribution extending up to the beam 
energy is observed with a 30% relative density. This tail 
must result from an instability occurring closer to the 
beam source. 

The low velocities of the electron holes (0.1-0.2ub ea m) 
and their lack of dependence on the beam velocity ex- 
cludes generation by a two-stream instability, which here 
predicts a velocity » 2^b e am/3. The two-stream in- 



stability generation mechanism was discarded for anal- 
ogous reasons by [H and [l8[. Similarly the electron 
holes or electrostatic whistler mode waves in this fre- 
quency range cannot be directly generated by a resonant 
(lu — fc||fbcam = 0) interaction with beam electrons or a 
fan instability based on anomalous Doppler-shifted res- 
onance (uj — k»vn = —LOce) with electrons from the tail. 
However the observed electron distribution carries a sub- 
stantial parallel current. Parallel currents in low- j3 plas- 
mas can drive unstable oblique electrostatic whistlers up 
to the frequency range of our measurements [l9| (for case 
b in this article). In our experiment a 62 eV 25% beam 
results in a mean electron velocity ks AAvtc, in which 
case this instability predicts a parallel phase velocity of 
2.2vTe, much closer to the measurements. Note however 
that the observed electron distribution has a background 
electron population not considered in [l9| . It also remains 
unclear if this instability can directly generate electron 
holes or they form later as the waves travel away from 
the beam center and the region of instability. 

In summary, for the first time a laboratory experiment 
has shown the generation of localized electrostatic struc- 
tures by the injection of a suprathermal electron beam. 
The scales and amplitudes of these structures are compa- 
rable to those derived from observations in various places 
in the magnetosphere, and their properties are consis- 
tent with electron holes. The electron holes are found 
on the high frequency (/ < f pe ) end of the electrostatic 
whistler mode which suggests that they have the same 
origin. We conjecture they result from a lower-hybrid 
instability driven by parallel currents. 
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